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The hydroxide anion plays an essential role in many chemical and biochemical reactions. But 
a molecular-scale description of its hydration state, and hence also its transport, in water is cur- 
rently controversial. The statistical mechanical quasi-chemical theory of solutions suggests that 
HO • [H2 0]3~ is the predominant species in the aqueous phase under standard conditions. This 
result is in close agreement with recent spectroscopic studies on hydroxide water clusters, and with 
the available thermodynamic hydration free energies. In contrast, a recent ab initio molecular dy- 
namics simulation has suggested that HO ■ [H20]4~ is the only dominant aqueous solution species. 
We apply adiabatic ab initio molecular dynamics simulations, and find good agreement with both 
the quasi-chemical theoretical predictions and experimental results. The present results suggest a 
picture that is simpler, more traditional, but with additional subtlety. These coordination struc- 
tures are labile but the tri-coordinate species is the prominent case. This conclusion is unaltered 
with changes in the electronic density functional. No evidence is found for rate-determining ac- 
tivated inter-conversion of a HO - [H20]4~ trap structure to HO - pHhOJa - , mediating hydroxide 
transport. The view of HO - diffusion as the hopping of a proton hole has substantial validity, the 
rate depending largely on the dynamic disorder of the water hydrogen-bond network. 



I. INTRODUCTION 

A pre-eminent challenge in liquid state physics is the 
understanding of aqueous phase chemical transforma- 
tions on a molecular scale. Water undergoes limited au- 
toprotolysis, which is enhanced in the presence of highly 
charged metals, such as Be 2+ 0, ■ Understanding the 
hydration and transport of the autoprotolysis products, 
H + and HO~, presents unique and interesting challenges 
for molecular-scale theories of solutions and for simula- 
tions. In this paper we focus on HO~(aq). 

Since H + and HO~ constitute the underlying aqueous 
matrix, it is not unreasonable to expect that their trans- 
port in water is different from the transport of other aque- 
ous ions. This anomalous diffusion of the H + (aq) and 
HO~(aq) has received extensive scrutiny over the years 
(for example, 0,0,0), but recently ab initio molecular 
dynamics (AIMD) capabilities have evolved to provide 
new information on the solution condition and transport 
of these species. Over a similar period of time, the sta- 
tistical mechanical theory of liquids (especially water) 
has also become usefully more sophisticated (for exam- 
ple, |6|]). These two approaches can be complementary, 
but in typical practice they remain imperfectly connected 
(but see BSBSElini) 

In an initial AIMD study [12, HCT(aq) was observed 
to be tetra-hydrated during the course of the ~6 ps long 
simulation. This complex had a lifetime of about 2-3 ps. 
An approximately square-planar configuration was noted 
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for this HO • [H20]4~ complex. That study hinted that 
transport occurred when HO ■ [H20]4~ converted to a 
tri-hydrated (HO • [H20]3~) species which has hydrogen- 
bonding arrangements similar to those in liquid water. 

A recent AIMD study 0] reinforced the notion 
of a HO • [H 2 0]4~ species exclusively dominating the 
equilibrium population distribution at infinite dilution. 
The proposed mechanism for HO~(aq) transport was: 
first the stable (and hence inactive) HO • [H20]4~ con- 
verts to the active HO • [H2 0]3~ species; then the H- 
bond between the anion and one of the ligating wa- 
ter molecules shortens, thus identifying a transient 
HO ■ [H20]~ species; a shared proton is transferred along 
that shortened bond, and a HO ■ [H20]3 _ species is re- 
constituted with the hydroxide identity now switched; 
this active HO • [PbOJs - species reverts back to an in- 
active HO • [H20]4~ species completing one transport 
event. Presumably, as in the lifetime of the 

HO • [H20]4~ species was 2-3 ps, but statistical charac- 
terization was sketchy. 

Discussions of a transport mechanism for HO _ (aq) 
typically focus on Agmon's 0, 0] extraction of an acti- 
vation energy for hydroxide transport from the temper- 
ature dependence of the experimental mobilities. Near 
room temperature that empirical parameter is about 
3 kcal/mol but increases by roughly a factor of two for 
slightly lower temperatures. As a mechanical barrier this 
value, about 5-6 krjT, may be low enough to require some 
subtlety of interpretation [16j ; the observed temperature 
sensitivity of the activation energy, and particularly its 
increase with decreasing temperature, supports that pos- 
sibility. We note that a standard inclusion of a tunneling 
correction would be expected to lead to a decrease of 
activation energy with decreasing temperature. 
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However, Ref. ^| framed the consideration of HO - 
transport in terms of classical transition state theory and 
extracted an activation energy from the gas-phase study 
of Novoa et al. [l7|. Ref. |13| also considered the im- 
portance of tunneling in lowering the barrier for proton 
transfer by performing path integral calculations. Their 
combined value of 3.1 kcal/mol was close to Agmon's es- 
timate |l4l Il5| . That experimental number does reflect 
the influence of the solution medium, and the Arrhenius 
plots are non-linear. Additionally, the earlier gas-phase 
studies 01 have shown that outer shell disposition of 
the fourth water molecule HO • [H2 0]3 - • H2O is lower 
in energy than HO • [H20]4^, in agreement with subse- 
quent experimental results [18j| : and, further, with a bar- 
rier of 2.5 kcal/mol for conversion HO • [H20]3 - • H2O to 
HO • [H20]4~, nearly 1.3 kcal/mol greater than the bar- 
rier for the reverse process. 

A follow-up classical AIMD study which treated 
a 1.5 M solution of KOD (1 molecule of OD - and 1 atom 
of K + in a 32 water molecule system) , estimated the dif- 
fusion coefficient of OD~ to be 2.1 A 2 /ps. A tunneling 
correction for the classical treatment would be expected 
to increase this rate. The experimental diffusion coeffi- 
cient for the light water OH - case is about 0.5 A 2 /ps at 
298 K @. 

The incongruities in |l2t Il3| were noted recently by 
H3 in AIMD simulations of deuterated NaOD and KOD 
hydroxide solutions with concentrations ranging between 
1.5 M to 15 M HJ. Interestingly, they observed that 
HO • [H2 0]3~ was well- represented in the population dis- 
tribution, and commented that their results |20| differed 
somewhat from those obtained in the previous ab initio 
study |l2Ul.'j| which reported HO • [H20]4~ to be the only 
dominant solvation structure. The distribution of hy- 
dration numbers was markedly influenced by the cations 
there [2pj. Those results do not permit a consistent ex- 
trapolation of HO~(aq) properties to infinite dilution, 
but HO ■ [H20]4 _ was just a prominent structure, not 
the only one. 

It is clear that earlier simulations have not resolved 
the most primitive question: What is the coordination 
state of HO~(aq) in water without extrinsic complica- 
tions? Speculations regarding the transport mechanism 
are somewhat premature in the absence of a clear un- 
derstanding of this coordination number question. Here 
we focus on that primitive question first, apply AIMD 
methods, and find that HO ■ [H20]3~ is a probable co- 
ordination number. We then discuss the agreement of 
the present simulation results with inferences based on 
(a) molecular theory pH l22l|. (b ) sp ectroscopic [Tsl ] and 
thermo-chemical measurements 1231 on hydroxide water 
clusters, (c) spectroscopic studies |24l l25l 1261 127| of aque- 
ous HO - , and (d) dielectric dispersion measurements of 
aqueous HO - [28j]. Discussion of the transport can be 
then framed within the dynamical disorder framework 

H^, HO, HH Q . 



II. RECENT THEORETICAL AND 
EXPERIMENTAL BACKGROUND 

Recent experimental and theoretical results have ad- 
dressed the issue of the coordination number of the aquo 
hydroxide ion. On the theoretical side, the statistical 
mechanical quasi-chemical theory of solution has been 
applied to HO - (aq) 0, |2^. This formally exact ap- 
proach, with roots in the work of Guggenheim [33L l34| 
and Bethe |35j , acquires approximations as applied. But 
for both hydration of ions in water and standard pack- 
ing problems, simple approximations have been proven 

effective 0, la (a S HE H3, HI, H3, E3, E3, 

In the quasi-chemical approach, the region around the 
solute is partitioned into inner and outer shell domains. 
The inner shell, where chemical effects are important, 
is treated quantum mechanically. The outer shell con- 
tributions can be assessed u sing classical force-fields or 
dielectric continuum models [lOj. The theory permits a 
variational check of the partition [E • Quasi-chemical 
studies jHll^ have firmly suggested that HO • [H 2 0] 3 - 
is the most prominent solution species 0, |2^ ■ Those 
results are insensitive to the choice of density functional 
or ab initio M0ller-Plesset 2 nd order perturbation tech- 
nique. 

Recent cluster spectroscopic studies ^8| have observed 
shell closure with formation of HO • [H20]3 - in the hy- 
dration of HO - . The fourth water molecule initiates an 
outer shell around this cluster. This identification of a 
shell-closure is in agreement with thermo-chemical mea- 
surements that show the same effect 23]. Thus theoret- 
ical considerations and experiments on hydroxide-water 
clusters concur on the significance of HO • [H2 0]3~ as the 
nominal full inner shell structure. 



III. AB INITIO MOLECULAR DYNAMICS 
SIMULATIONS 

The AIMD simulations were carried out with the 
VASP |44l |45| simulation program wherein the Kohn- 
Sham equations are solved by usual matrix methods. 
Thus this corresponds to adiabatic dynamics, with tol- 
erance set by the convergence criterion for the electronic 
structure calculation. The system comprises a hydrox- 
ide anion in a periodic cube of 32 water molecules. The 
box size was 9.8788 A consistent with the experimental 
partial specific volume of the HO~(aq) 46]. This sys- 
tem was initially thermalized by about 10 ps of classical 
molecular dynamics (using SPC/E [4jJ potentials) with 
a temperature of 300 K using velocity scaling. The dom- 
inant HO - coordination number during this phase was 
n = 5. 

In the first AIMD simulation, RUN1, we used a gen- 
eralized gradient approximation, PW91 HE, EH' to elec- 
tron density functional theory. Ultrasoft pseudopoten- 

tiais (us-pp) mim 

for oxygen were used to describe 
the core-valence interaction, and an ultrasoft pseudopo- 
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tential was used for the hydrogen atoms as well|44l |45|. 
The valence orbitals were expanded in plane waves with 
a kinetic energy cutoff of 395.76 eV. RUN1 used a 1 fs 
timestep. SCF convergence was accepted when the en- 
ergy difference between successive iterations fell below 
1 x 1CP 4 eV. Initially 1.5 ps of AIMD simulations were 
performed at 300 K using velocity scaling to maintain 
the stated temperature. The coordination number at the 
end of this 1.5 ps simulation was n — 4. The system was 
then equilibrated for 3.2 ps in the NVE ensemble. Then 
random velocities were reassigned to give a temperature 
of 300 K, and statistics collected for another 8.2 ps in 
the NVE ensemble. The mean t empe rature was 332 K. 
The relative energy fluctuation, V SE 2 /\E\, was 8.4- 10~ 5 . 



The drift in the relative energy was about 8-10 6 
These values appear quite reasonable |52l |. 



ps" 



To evaluate the role of chosen density functionals and 
the timestep, we took the terminal configuration from 
RUN1 and replaced all the hydrogen atoms by deu- 
terium, thus simulating the classical statistical mechan- 
ics of aqueous DO~ in D 2 0. The timestep was also 
cut in half to 0.5 ps. After a short NVE run with 
the PW91 density functional and US-PP pseudopoten- 
tials, we started RUN2 and RUN3. In both these runs 
the pseudopotential treatment of atoms was replaced 
by the projector augmented-wave [53L 154^ (PAW) treat- 
ment, which is thought to handle difficult cases involving 
large electronegativity differences with "exceptional pre- 
cision" |54j. Further, in molecular bonding problems, this 
method is about as accurate as local basis (such as Gaus- 
sian orbital) methods H^. In RUN2 and RUN3, the SCF 
convergence was accepted when the energy difference be- 
tween successive iterations fell below 1 x 10 -6 eV. (For 
comparison, this is an order of magnitude smaller than 
the "tight" convergence in the Gaussian |5y| suite of pro- 
grams.) 




FIG. 1: R, is the radius of the observation volume centered on 
the hydroxyl oxygen. 9 and <j> identify the angles that specify 
the directionality of the H-bond to water. 
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FIG. 2: RUN1 coordination number and temperature versus 
time. R = 2.5 A. The block averaged temperature is shown 
with the solid line. The mean temperature is 332±22 K. The 
short vertical bars at the n = 3.5 level flag hydrogen exchange 
events, which also changes the identity of the hydroxyl. Note 
that many hydrogen exchange events occur without interces- 
sion of the n=4 configuration. 



IV. DISCUSSION OF AIMD RESULTS 



RUN2 employed the PBE functional H3 which is sim- 
ilar to the PW91. RUN2 is not considered further here 
because it produced a predominant coordination number 
of n = 3, just as RUN1 did. RUN3 employed the revised 
PBE functional j5^. The system was equilibrated for 
5.9 ps and a further production run of 5.9 ps was con- 
ducted. The mean temperature was 313 K. VSE 2 /\E\, 
was 2.0 • 10 -5 . The drift in the relative energy was about 
5 • 10- 6 ps^ 1 . 

FIG. n introduces the geometric notation used in ana- 
lyzing the coordination of HO" (aq). FIG. |U and FIG. 0| 
show the coordination number at each time for RUN1 
and RUN3, and also the instantaneous temperature ob- 
served. Radial distribution functions for those two cases 
are shown in FIG. Hand FIG. H TABLE [Q presents the 
average fractional coordination number populations for 
RUN1 and RUN3, and TABLE HU records averages of H- 
bonding angles using the notation of FIG.[T] 



For RUN1, it is clear that n — 3 is the predominant 
state (Figs. Hand Eland TABLE 0). The radial distribu- 
tion functions shown in FIG. (bottom panel) establish 
that R < 2.5 A is a reasonable, inclusive, even permis- 
sive selection criterion; see also |59| . TABLES H] and |n] 
provide guidance on whether many of these n = 4 con- 
figurations should be excluded as not hydrogen-bonded. 
Note that mediating n = 3 to n = 4 inter-conversions are 
not required for exchange of the identity of the hydroxide 
oxygen. RUN3 below provided the same observation. 

For RUN3, using the R < 2.5 A criterion, we find 
about equal populations of n = 3 and n — 4. Tighten- 
ing this criteria by 0.25 A drops n—A population by 40% 
(TABLE HJl relative to n=3. It is apparent from TA- 
BLE [H] that many of these n = 4 states are not square- 
planar. TABLE [I] shows that even a permissive 9 > 80 
cutoff excludes many of those n = 4 cases. The configu- 
rations thereby excluded are on the 'forward' side of the 
hydroxide- water complex, above the plane defined by the 
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FIG. 3: Density distribution of water oxygen and proton 
around the hydroxyl oxygen for RUN1 which utilized PW91. 
The distributions of the neighboring atoms are also separated 
into contributions according to distance-order. The hydrogen 
of the nominal HO chemical bond, otherwise the nearest H, 
isn't included in this distance-ordering. 
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FIG. 5: Density distribution of water oxygen and proton 
around the hydroxyl oxygen for RUN3 which utilized rPBE; 
otherwise as in Fig. El 



TABLE I: Relative populations Xi — Xi /x 3 for RUN1 and 
RUN3 and different selection (FIG. criteria. The quasi- 
chemical theory results are denoted by QC. 



Method 


Criteria xi £2 £3 xa 


QC 


0.03 0.26 1.0 0.0 


PW91/RUN1 
PW91/RUN1 


R < 2.5 1.0 0.14 
R < 2.5,6* > 80,0 > 150 -- 0.03 1.0 0.08 


rPBE/RUN3 
rPBE/RUN3 
rPBE/RUN3 
rPBE/RUN3 


R < 2.5 1.0 1.02 
R < 2.25 - 0.01 1.0 0.60 
R < 2.5,6* > 80 - 0.01 1.0 0.43 
R < 2.5,9 > 80,0 > 150 -- 0.07 1.0 0.36 



FIG. 4: RUN3 coordination number and temperature ver- 
sus time. The mean temperature is 313±21 K. Rest as in 
FIG. [5] The dashed line applies to the selection criterion in- 
volving R < 2.5 , 8 > 80 , (f) > 150. Note that many hydrogen 
exchange events occur without intercession of the n=4 con- 
figuration. 



hydroxyl oxygen and the shortest OH bond 

The radial distributions (FIG. El and FIG. 0) decom- 



posed according according to the distance-order of atoms 
surrounding the hydroxide oxygen, denoted O*, are sim- 
ilarly interesting. Though the statistical quality is mea- 
ger, the conventional O atom second shell begins at a 
distance roughly \[2 times the radial location of the first 
shell. The fourth-nearest oxygen atom, qualitatively de- 
scribed, builds a shoulder on the outside of the principal 
maximum of the 0*0 radial distribution functions. Note 
that the contributions from the nearest three protons and 
the nearest three oxygen are concentrated, and that those 
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TABLE II: Mean angles 9 and cj>, in degrees, defined in Fig.0 
for RUN1 and RUN3. The four values pertain to the four 
nearest coordinating protons ordered in distance from the an- 
ionic oxygen. The values of (9j), j=l,2, and 3 are consistent 
with classic tetrahedral geometry (109.5°). These values are 
in good agreement with angles (110°) obtained from the op- 
timized HO • [H20]3~ cluster. {64) is different; 64 is typically 
located closer to the equatorial plane, but with bigger statisti- 
cal dispersion. (64) is also different from the angle (116°) ob- 
tained from the optimized HO • [H20]4~ cluster. The angles 
((j>j) indicate that the coordinating OH bond is not collinear 
with the 0*0 vector, and this is consistent with the cluster 
results. Note specifically, that the water oxygen atom deter- 
mining the angle cj>j doesn't corrrespond uniquely to a partic- 
ular distance order for oxygen atoms; this angle is defined by 
the distance-ordering of the hydrogen atoms, and the oxygen 
atoms to which those hydrogens are directly bonded. 



(fli) (O2) (fe) (64) 

PW91/RUN1 107.9±9.2 108.1±10.4 107.0±11.8 97.4±32.8 

rPBE/RUN3 109.1±9.9 108.0±11.9 103.4±13.1 91.8±24.5 

(<M (<M (<fe) (<M 

PW91/RUN1 169.2±9.9 168.1±6.3 166.7±7.1 84.2±38.8 

rPBE/RUN3 169.8+5.4 168.0+6.0 165.0+7.9 124.9+45.8 



protons are about lA nearer the hydroxide oxygen. In 
contrast, the contributions from the fourth-nearest pro- 
ton and fourth-nearest oxygen are diffuse and overlap- 
ping; the contribution from the fourth-nearest proton 
is not always inside the contribution from the fourth- 
nearest oxygen. (Of course, those atoms needn't be di- 
rectly bonded.) These observations suggests again that 
the fourth-nearest water molecule is not always partic- 
ipating in a conventional, specific H-bond but is often 
non-specifically arranged. This description as a whole is 
consistent with the fact that the tri-coordinated species is 
a prominent species though the tetra-coordinated species 
is present to some extent also. 

The (Qj) results of TABLElLTldocument the interesting 
point that three nearest coordinating protons are physi- 
cally equivalent and approximately disposed towards the 
corners of a tetrahedron. The fourth-nearest proton is 
distributed broadly about the plane containing the hy- 
droxyl oxygen and perpendicular to the OH chemical 
bond. 

Figs. El and show that the nearest water-oxygen is al- 
ways near 2.45+0. 1(2ct) A of the hydroxyl oxygen. This 
distance is close to the 0-0 separation in the calculated 
gas-phase structure of HO • [H2O] - , 2.46A. We conclude 
that HO - [H20]~ is a prominent sub-grouping in the 
HO • [H2 0] n ~ (n — 2,3,4) species. The recent theoreti- 
cal and experimental work noted in Sec. HT1 suggests that 
the structure HO • [H20]~ does not provide as natural a 
description of the thermodynamic hydration free energy 
as does the species HO ■ [F^O]^ that follow from a more 
conservatively defined inner-shell. It seems likely that an 



inner-shell definition designed |4£j to better isolate this 
HO • [H20]~ sub-grouping would be excessively compli- 
cated. This is one reason why xi is not unambiguously 
separated by the results of TABLE [I] 

The identification of HO • [H20]~ as a prominent sub- 
grouping agrees with spectroscopic studies on concen- 
trated hydroxide solutions. The IR and Raman spectra of 
concentrated hydroxide solutions have been interpreted 
in terms of HO • [H 2 0]~ as a p rincipal structural possi- 
bility for those systems |13, |25| 0, H3 • But the present 
observations suggest that it is the principal participant 
in the proton conductivity. 

This HO • [F^O] - sub-grouping also concisely resolves 
the high "effective" (not microscopic) hydration num- 
bers extracted from dielectric dispersion measurements 
|28j. A 'super'-grouping of hydrated HO - [H 2 0]~, one 
involving several more water molecules, could well be 
relevant to the time scale of the measurement, a pos- 
sibility also suggested by Agm on [lp. Then a dominat- 
ing HO ■ [H 2 0]4~ species |iH, liH EfJj is not a conclusion 
necessary to the resolution of experimentally obtained 
"effective" hydration numbers. 

The diffusion coefficient of HO - is calculated to be 
3.1 A 2 /ps at 332 K using PW91 for HO", and 1.1 A 2 /ps 
at 313 K using rPBE for DO~. For comparison, Zhu and 
Tuckerman's results for DO~ gives 2.1 A 2 /ps presumably 
at 300 K [l^. The experimental value, calculated from 
mobility data, is about 0.5 A 2 /ps for HO - 0. Perhaps 
the agreement with experiment of this rPBE result is ac- 
ceptable, but all these values agree only roughly with ex- 
periments. An important background point is that eval- 
uation of diffusion coefficients typically requires longer 
simulation times, not readily accessible here by AIMD 
simulation. 

Rationalizations why the simulation rates arc higher 
than observed experimental rates would be highly spec- 
ulative. A quantum mechanical treatment of the water 
matrix, rather than a classical one, would imply less or- 
der, and perhaps that would lead to slower rates. Or it 
might be that the electron density functionals used here 
lead to excessive prominence of the HO • [H 2 0]~ sub- 
grouping, and that leads to rates that are too high by 
comparison with experiment. The fact that the empirical 
activation energies increase with decreasing temperature 
remains an unexplained point, apparently of qualitative 
significance. 

The transport observed here involved the movement of 
a proton-hole between two tri-hydrated species, which is 
also consistent with the identification of the HO ■ [H20]~ 
sub-grouping and the IR spectra. Second-shell rearrange- 
ments do occur, but all-or-nothing breaking and reform- 
ing of a hydrogen bond is not necessary |l5|. There will 
certainly be rearrangements as the hole settles into its 
new place. The hole-hopping proposal for HO~ trans- 
port, as discussed by Bernal and Fowler Q, Hiickel 0, 
and later Stillinger [5j, has substantial validity. 

We note that AIMD studies on pure liquid water 0] 
under conventional thermodynamic conditions show that 
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PW91 and PBE predict more strongly structured liquid 
water compared to experiment, and rPBE softens the 
structure of liquid water simulated on that basis. The 
excess chemical potential evaluated in that study |43j . 
on the basis of AIMD results and the quasi-chemical the- 
ory, indicates that PW91 binds liquid water too strongly, 
whereas rPBE softens the binding. The excess chemical 
potential of water at 314 K using the rPBE functional 
was —5.1 kcal/mole (compared with —6.1 kcal/mole ex- 
perimentally), whereas PW91 gives —12 kcal/mole at 
330 K. Despite these differences, the hydration struc- 
ture of HO _ (aq) and the qualitative transport pattern 
are similar in these two cases. Therefore, the above com- 
parisons can be more optimistic for providing the correct 
direction to frame discussions on HO~ transport. 

V. CONCLUSION 

Three distinct lines of investigation, theory [221 . ex- 
periments 0, |2^| , and the present simulations |2l|. con- 
verge on the common view that HO • [H2 0]3~ is a promi- 
nent, likely even dominating coordination structure for 



HO~ (aq); this is the most primitive issue underlying 
current speculations regarding HO~ in aqueous solu- 
tions. The present simulation results suggest, in addi- 
tion, that the coordination number distribution is la- 
bile, includes less specifically structured n=A possibili- 
ties, and that HO • [H20]~ is a prominent sub-grouping 
within larger inner shell structures. This latter point is 
consistent with in terp retations of aqueous phase spectro- 
scopic |2J, |23, |26|, I23 results, and also with proposals of 
high-effective solvation numbers on the basis of dielectric 
dispersion measurements pgj . 
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